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ABSTRACT 

We examine, in a general manner, the role played by the initial density and 
velocity distributions in the gravitational collapse of a spherically symmetric inho- 
mogeneous dust cloud. Such a collapse is described by the Tolman-Bondi metric 
which has two free functions: the 'mass-function' and the 'energy function', which 
are determined by the initial density and velocity profiles of the cloud. The collapse 
can end in either a black-hole or a naked singularity, depending on the initial pa- 
rameters characterizing these profiles. In the marginally bound case, we find that 
the collapse ends in a naked singularity if the leading non-vanishing derivative of 
the density at the center is either the first one or the second one. If the first two 
derivatives are zero, and the third derivative non-zero, the singularity could either 
be naked or covered, depending on a quantity determined by the third derivative 
and the central density. If the first three derivatives are zero, the collapse ends 
in a black hole. In particular, the classic result of Oppenheimer and Snyder, that 
homogeneous dust collapse leads to a black hole, is recovered as a special case. 
Analogous results are found when the cloud is not marginally bound, and also for 
the case of a cloud starting from rest. A condition on the initial density profile is 
given for the singularity to be globally naked. We also show how the strength of the 
naked singularity depends on the density and velocity distribution. Our analysis 
generalizes and simplifies the earlier work of Christodoulou and Newman [4,5] by 
dropping the assumption of evenness of density functions. It turns out that relaxing 
this assumption allows for a smooth transition from the naked singularity phase to 
the black-hole phase, and also allows for the occurrence of strong curvature naked 
singularities. 
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I. Introduction 

In the framework of general theory of relativity, the end state of gravitational 
collapse of a sufficiently massive body is a gravitational singularity. However, it is 
not known whether such a singularity will be naked or covered by an event horizon. 
The well known Hawking-Penrose singularity theorems provide no information on 
this issue. The cosmic censorship hypothesis essentially states that in the general 
theory of relativity, the end state of gravitational collapse is always a black hole: the 
gravitational singularity must necessarily be covered by an event horizon. Whereas 
no proof for this hypothesis is known, many counter-examples have been studied 
recently. It is clearly important to investigate specific dynamical collapse scenarios 
within the framework of gravitation theory for the final state of collapse. Such a 
study will be of great help in establishing whether or not this hypothesis is correct. 

An important class of collapse scenarios available in this connection is the 
dynamical evolution of a spherically symmetric cloud of pressureless dust. Dust 
collapse has been studied by many workers, perhaps the first being the work of Op- 
penheimer and Snyder [1] , which demonstrated that the collapse of a homogeneous 
spherical dust object ends in a black-hole. This leads to a view, reflected by the 
cosmic censorship hypothesis, that spherical collapse will always end in a black-hole. 
However, such a view does not seem to be supported by a detailed analysis of spher- 
ical collapse. Departures from the Oppenheimer-Snyder model come in the form 
of introducing inhomogeneities in the initial distribution of matter [2], and also in 
the form of changing the equation of state. It has been shown by various authors 
[3-9] that the introduction of inhomogeneities in the matter distribution can give 
rise to the occurrence of a naked singularity at the center in spherical gravitational 
collapse. Also, equations of state, different from dust, have been investigated in this 
context [10-15] (see e.g. Ref.[14] for a review on these and related developments). 
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One may try to take the viewpoint that the censorship hypothesis has to be suit- 
ably modified to exclude such examples of naked singularities. Such an approach 
requires great care because in some of these studies it is seen that depending on 
the initial conditions, the collapse ends in a naked singularity or a black hole, and 
there is a smooth transition from one phase to the other [16]. 

Given that general relativity does seem to predict the occurrence of naked 
singularities in gravitational collapse, it would be important to figure out the initial 
conditions which give rise to same. In particular, one would like to find out the role 
played by the initial distribution of matter. For instance, in the case of dust, which 
density distributions lead to a naked singularity, and which ones to a black hole? 
Could it be that the initial distributions which lead to naked singularities are not 
astrophysically realizable? 

We have started a program of this nature, and in this paper we consider the 
effect of various initial density and velocity distributions on the final fate of gravita- 
tional collapse of spherically symmetric inhomogeneous dust. The principal reason 
for considering dust, even though it does not take into account the pressures which 
might be important in the later stages of collapse, is to pave the way for consid- 
eration of more general forms of matter. The occurrence of naked singularities in 
spherical collapse, when a general form of matter is considered has been demon- 
strated [15], although the explicit role of the initial matter distribution remains to 
be investigated. 

It is relevant to note here that the role of initial density distribution in dust 
collapse was earlier considered by Christodoulou [4] and Newman [5]; however, their 
discussion was restricted by the assumption of evenness of density and metric func- 
tions. From an astrophysical viewpoint, we merely assume the analyticity of all 
the physical functions concerned, but drop the restrictive assumption of evenness. 
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An interesting structure in the transition from the naked singularity phase to the 
black hole phase then emerges when the assumption of evenness is relaxed. The 
present work simplifies and adds clarity to the earlier discussions of Christodoulou 
and Newman, apart from generalizing it. For this purpose, we use here the method- 
ology developed in [9] for the study of collapse, which examines the nonspacelike 
geodesies of the spacetime near the singularity. One constructs there an algebraic 
equation; the existence of the real positive roots for the same characterizes the for- 
mation of a naked singularity or black hole in the collapse. The solutions which 
correspond to the existence of future directed null geodesies terminating in the past 
at the singularity could then be determined. When such solutions exist, the sin- 
gularity is naked. This yields constraints on the physical parameters describing 
the initial data for the collapse, if the singularity is to be naked. It is straightfor- 
ward then to translate these restrictions into constraints on the density and velocity 
distribution of the dust cloud, as shown here. 

Here is a summary of the main results in this paper. The gravitational collapse 
of spherically symmetric inhomogeneous dust is described by the Tolman-Bondi 
models [2]. This metric has two free spatial functions: the 'mass function' F(r) and 
the 'energy function' f(r). These functions get determined once the initial density 
and velocity profiles of the collapsing dust cloud are given. We find that the nature 
of the singularity is determined by the behavior of the initial density and velocity 
profiles near the center. In the marginally bound case, (i.e. / = 0), it is found 
that the collapse leads to a locally naked singularity if the leading non-vanishing 
derivative of the density at the center is either the first or the second derivative. 
Also, the singularity is gravitationally weak. If the first two derivatives are zero 
but the third one is non-zero, the collapse ends in either a black hole or a naked 
singularity, depending on the value of a quantity determined by the third derivative 
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and the central density. In this case we get a strong curvature singularity. If the first 
three derivatives are zero, the collapse ends in a black hole. Analogous results are 
found when the collapse is not marginally bound, and in particular, when the cloud 
starts its collapse from rest. In this manner, we unify the results of Christodoulou 
[4], Newman [5] and Dwivedi and Joshi [8] on spherical dust collapse. 

The plan of the paper is as follows. In Section II we recall the equation which 
determines whether the singularity resulting from collapse is naked or covered. In 
Section III we use this equation to determine the final fate of a collapsing cloud. 
The cases of a marginally bound cloud, a cloud starting from rest, and a cloud 
which is not marginally bound are considered separately. In Section IV we examine 
issues relating to the strength of the singularity, and its global nakedness; this is 
followed by a general discussion. 

II. Naked singularity in the Tolman-Bondi gravitational collapse 

In this section we set the terminology and also recall, for the sake of complete- 
ness, the equation derived in [9], which determines whether a naked singularity can 
form in the gravitational collapse of a spherically symmetric dust cloud. 

A spherically symmetric inhomogeneous dust cloud is described by the Tolman- 
Bondi metric, which in comoving coordinates is given by 

ds 2 = -dt 2 + -^rj dr2 + R2 ( d02 + sin2 9d ^- t 1 ) 

/ is an arbitrary function of the comoving coordinate r, satisfying / > —1. R(t, r) is 
the physical radius at time t of the shell labeled by r, in the sense that 4nR 2 (t, r) is 
the proper area of the shell at time t. Prime denotes partial derivative with respect 
to r. 

The energy-momentum tensor is that of dust given by T iJ = e5 l t 5 J tl where e(t, r) 
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is the energy density of the cloud, and the Einstein equations for this metric are 



R 2 = ^ + /■ (3) 

(We have set 87rG/c 4 = 1). F is an arbitrary function of r; dot denotes partial 
derivative with respect to time. Since we are concerned with collapse, we assume 
R < 0; the equality sign holds at the initial epoch for a cloud whose collapse starts 
from rest. Eqn.(3) can be integrated to obtain R as a function of t and r, given 
implicitly by the relation 

t - toW = - fl3 % /fl/F) , (4) 
with G(y) a positive function having the range 1 > y > — oo and given by 



a^_vS l>y>0, 



y 3/2 y 



G(y) = <f _ y = 0, (5) 

Here to( r ) is a constant of integration which we determine by noting that there is 
a scaling freedom in the choice of r. Using this scaling freedom, we require that at 
the starting epoch of collapse, t = 0, we have i?(0, r) = r. From Eqn. (4) it then 
follows that 

to(r) gives the time at which the physical radius of the shell labeled by r becomes 
zero, and hence the shell becomes singular. Unlike the case of collapse of a homo- 
geneous dust cloud, in the inhomogeneous case different shells become singular at 
different times. 



If there are future directed radial null geodesies coming out of the singularity, 
with a well defined tangent at the singularity, then the quantity R' must tend to 
a finite limit in the limit of approach to the singularity in the past along these 
trajectories. This follows from (1). The partial derivative R' will appear in these 
calculations, and it is necessary to write it in such a way that it approaches a well- 
defined limit as the singularity is approached. Using Eqns. (3)-(6) we can write it 
in the form 



R' = r a ~ x 

„a-l 



(r, - 0)X + [0-(v- l(3)X^G(-PX)] [P + i] 1 ^ 



(7) 



r a - L H(X, r) 



where 



X = (8a) 
tF' 

V = W) = -pr, (86) 



/ 

p = p(r) = r -l, (8d) 



P = pr a ~\ (8e) 

A=^ (8/) 
t^/X 1 + 0-r, (y-lP)G(-p) 

^ — r a-l _|_^l/2 r 3(a-l)/2 + r 3(a-l)/2 ) • 

The function (3{r) is defined to be zero when / is zero. The factor r Q has been 
introduced here for the sake of convenience in examining the structure of the naked 
singularity. The exact value of the positive constant a > 1 is to be determined and 
will depend on the specific model considered, as we show in the following sections. 

The points (t , r) where a singularity R(t , r) = occurs are given by Eqn. (4). 
This corresponds to the physical situation when the matter shells are crushed to 



zero radius. The singularity corresponding to r = is called the central singularity. 
The singularity is naked if there are future directed non-spacelike curves in the 
space-time with their past end-point at the singularity. Here we restrict ourselves 
to the study of future directed radial null geodesies. In order to check whether the 
singularity can be naked, we examine the null geodesic equations for the tangent 
vectors K a = dx a /dk, where k is an affine parameter along the geodesies. For the 
radial null geodesies these are 



dt V 



K 1 = — = 



dk R' 



Rr= dr = Vy/T+7 
dk RR' 

The function V(t, r) satisfies the differential equation 



(9) 
(10) 



hP 

dk 



VT+fR' _ R_ 
RR' R? 



R 



= 0. 



Writing the concerned quantities in terms of the variables R and u, where u = r a , 
the geodesic equation is written as: 

1 



dR 

du 



ar 



a-l 



■ dt , 
R— + R 

dr 



1 - 



(11) 



a 



VT+7 

If the outgoing null geodesies are to terminate in the past at the central singularity 
at r = 0, which occurs at time t = to at which R(t , 0) = 0, then along these 
geodesies we have R — > as r — > 0. In terms of variables u and R, the point u = 0, 
R = is a singularity of the above first order differential equation. For an outgoing 
null geodesic dR/du must be positive. 

As explained in [9] , the nature of the singularity is to be understood by exam- 
ining the behavior of the characteristic curves of the differential equation (11) in 
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the vicinity of the singularity. If these characteristics terminate at the singularity in 
the past with a definite tangent, this behavior is determined by the limiting value 
of the quantity X, which is defined by X = R/r a = R/u at R = 0, u = 0. If the 
null geodesies meet the singularity with a definite value of the tangent, then using 
(11) and l'Hospital rule we get, for the value of Xq, 

Xq = lim — = lim — — 

R^0,u^0 u R^0,u^0 du 

= lim U(X,u) = U(X ,0). 

R— >0,u— >0 

If a real and positive value of X satisfies the above equation the singularity will 
be naked. If no such real positive roots exist then the singularity cannot be naked, 
and the collapse ends in a black hole. 
We write Eqn. (12) as 

V(X ) = (13) 



where 

\ Vf + A /X 
v'T~7n 



V(X) = U(X, 0) -X = 



H{X '°'>-X. (14) 



a 



We have introduced the notation that a subscript zero on any function of r denotes 
its value at r = 0. 

The constant a represents the behavior of singular geodesies near the singu- 
larity. In fact we can write R = Xor a in the neighborhood of the singularity, Xq 
being the real and positive root of Eqn. (13). We also make the assumption that 
Xq is finite. This is reasonable because if it were infinite, a suitable redefinition will 
again make it finite. The constant a is determined uniquely by the requirement 
that 0(r) does not vanish or go to infinity as r — > in the limit of approach to the 
singularity along any X = constant direction. This is the defining condition for a. 
This condition ensures within our framework that the quantity H(X, 0) will not be 
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identically zero or infinite in Eqns. (11) and (14), and hence these equations will 
have a meaningful interpretation. (The only case in which O is identically zero is 
that of homogeneous collapse, for which 77 = 3 and f3 = 2 or 0.) 

III. The end state of dust collapse 

We shall now characterize in this section the formation or otherwise of a naked 
singularity in gravitational collapse in terms of the given initial density and velocity 
profiles for a dust cloud. The initial state of the spherically symmetric dust cloud 
is described in terms of its density and velocity profiles specified at an initial epoch 
of time from which the collapse commences. The density and velocity distributions 
together determine the free functions F(r) and f(r) as follows. We denote by p(r) 
the density distribution e(0, r) of the cloud at the starting epoch of collapse, t = 0. 
Since we are using the scaling R = r at t = 0, this is the true physical density of 
the cloud. From Eqn. (2) it then follows that 

F(r) = J p(r)r 2 dr. (15) 

Once F(r) is determined, we use it in (3), along with the initial velocity profile, 
v 2 (r) = -R 2 (0, r), to determine the function f(r). 

The question of interest is: given an initial density and velocity distribution 
what will be the end state of the collapse? In order to answer this question we 
assume that the density p(r) can be expanded in a power series about the central 
density po : 

00 n 

P^) = Y.Pn~v ( 16 ) 

n=0 

where p n stands for the nth derivative of p at r = 0. This gives, using (15), that 

00 

F(r) = Y,F n r n+3 , (17) 

n=0 
11 



where F n = p n j (n + 3)n\. 

A clarification should be made here on the nature of these power series expan- 
sions. Newman [5] and Christodoulou [4] assumed in their analysis that the metric 
and density, on a surface of constant time, considered as functions of r, are extended 
analytic. That is, after extending the coordinate r to include the range (— oo,0), 
they assumed these functions to be analytic on the whole real line. Further, they 
assumed these functions to be even, which automatically ensures that if the func- 
tion is analytic in the range zero to infinity, it will be extended analytic. Here, 
we demand that these functions be analytic only in the range zero to infinity, and 
we do not assume them to be even. Dropping the evenness assumption introduces 
important new qualitative features in the collapse scenario. It should be evident 
that there is no physical requirement as such that these functions be even. 

Given the density profile, we now determine a and Go as these quantities enter 
the equation for the roots, described in the previous section. For this, we need the 
expansion for the function r\ = rF' /F: 



We are interested in the behavior of r\ near the center. Hence, by using a binomial 
expansion for the denominator we are able to write r\ as follows: 



V = 



En = o(" + 3)F w r' 

p r n+3 



(18) 



r](r) = 3 + r]ir + r] 2 r 2 + ^r 3 + (9(r 4 ). 



(19) 



The leading order term is 3 and the next three terms are given by 



F 1 2F 2 Fl 



T t 1 = ~TP~ 




■2 ' 



3^3 Ff 3FiF 2 



(20) 



m = 
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It turns out that we will need the explicit form of terms up to order r 3 only. If all 
the derivatives p n of the density vanish, for n < (q — 1), and the qth derivative is 
the first non- vanishing derivative, then Tf], the qth term in the expansion for 77, is 

T* = ^*r«. (21) 

Here q takes the value 1, 2 or 3. (T° = 3). In this case, we can write rj(r) as 

v ( r ) = 3 + ^r q + 0{r q+1 ). (22) 

We first consider the marginally bound case / = 0, which gives that v 2 (r) = 
F(r)/r. Thus / = automatically implies that an initial velocity has been imparted 
to the cloud. The case of a cloud that is not marginally bound is considered next. 
This is followed by an examination of the cloud which starts from rest. 

A. Collapse of a marginally bound cloud 

From Eqn. (8), we take the expression for 0, with (3 = p = 0, and use 77 from 
(22) to get 

e - - qF i r<1 (23) 

3r 3(a-l)/2 3Fq r 3(a-l)/2' 

correct to order q in 77. The constant a is to be determined by the requirement that 
60, the limiting value of 6 as r — > 0, should not be zero or infinite. Clearly, for this 
determination, the leading non-vanishing term is the relevant one, and we get 

a = 1+ l- e °=-i- < 24 > 

All the higher order terms in G will go to zero as r — > 0, for this choice of a. 

Next, we return to Eqns. (13) and (14) for the roots. For the present case 
/ = 0, and using t]q = 3, these equations reduce to 

nx„)=0, V(X) = l{l-yf}{* + °L}-X, (25) 
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with 6 given by Eqn. (24). The limiting value of the function A = F/r a is found 
using (17) and (24) to be 

rO, q<3 

Ao = < F , q = S- (26) 
I oo, q > 3 

Since Ao takes different values for different choices of q, the nature of the roots 
depends on the answer to the following: which is the first non-vanishing derivative 
of the density at the center? Hence we analyze the various cases one by one, below. 

(i) The first derivative of density at the center, p\, is non-zero: 

In this case, q = 1, a = 5/3, and Eqn. (25) gives 

Y 3/2 _ F 1 3 pi 

We assume the density to be decreasing outwards, so p\ < 0, hence X will be 
positive and the singularity is naked. 

(ii) pi is zero, p2 is non-zero: 

In this case we have g = 2, a = 7/3, and from Eqn. (25) that 

Xl /2 = ~^L = -^. (28) 
2F 20p 1 ' 

Once again, we take p2 < 0, which means that X is positive and the singularity is 

naked. 

(iii) Pl = p 2 = 0, p 3 ^ 0: 

This perhaps is the most interesting case. We have q = 3, a = 3, A = F , and 
©o = —F3/F0 = —ps/12po. We take ps < 0. Since Ao 7^ 0, the equation (25) to be 
examined for real positive roots has a non-trivial structure. Substituting X = F x 2 
in (25) gives 

2x 4 + x 3 + £z-£ = 0, (29) 

where we have put £ = F 3 /F^ 2 = \^3p s /4p^ 2 . Using standard methods it can 
now be shown that this quartic has positive real roots if and only if £ < £2 or 
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£ > £i , where £1 = -.0096 and £ 2 = -25.9904. For £ 2 < £ < £i all the roots are 
imaginary. Thus, it follows that the singularity will be naked if and only if £ < £ 2 or 
£ > £i- When this condition is not satisfied the collapse ends in a black-hole. This 
is precisely the case considered by us in our paper [16], where a phase-transition 
from the black-hole phase to the naked singularity phase was found. The quartic 
in that paper looks different because a different scaling relating R and r was used. 
However it is easy to show that the parameter a defined in that paper (no relation 
with the a here) is related to the £ here by a = — 1/6£, so that when the condition 
on £ is translated into a condition on a, the naked singularity arises when a lies 
in the range 6.4126 x 10~ 3 < a < 17.3269. This is the same condition as derived 
in [16]. This recovers these earlier results here while working in a different scaling. 
The class of density profiles considered in [16] is indeed one where the first and 
second derivatives of the density are zero at the center, while the third derivative 
is non-zero, 
(iv) pi = p 2 = p 3 = 

When the first three derivatives of the density are zero at the center, then q > 4, 
a > 11/3, and Ao is infinite. It is easy to check in this case that a positive value of 
Xq cannot satisfy equation (25) and the collapse ends in a black-hole. In particular, 
the homogeneous collapse is seen to lead to a black-hole, as all the derivatives of 
the density are zero everywhere inside the cloud, including at the center. Thus we 
recover the results of Oppenheimer and Snyder [1] as a special case. 

The role of inhomogeneity is thus brought out clearly by our considerations 
here. If the leading non-vanishing derivative of the density at the center is either 
the first or the second one, the singularity is naked. If the leading non-vanishing 
derivative is the third one, the singularity is naked or censored, depending on the 
actual value of the ratio £ = ^/3ps/4p^ 2 . If the first three derivatives are zero, and 
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the first non-vanishing derivative is fourth or higher, the singularity is covered and 
the collapse ends in a black-hole. 

B. Collapse of a cloud that is not marginally bound 

We now include the energy function / in the analysis, which allows, in particu- 
lar, the collapse of a cloud starting from rest. The expansion for the functions p, F 
and r\ is the same as given above. The new functions that appear in the equation 
(13) for the roots are (3 = rf'/f, p = rf/F and P = pr a ~ l and we must evaluate 
their limits as r — > 0. Moreover, the limit 6 will now be different, because it will 
depend on /. 

As mentioned before, the function f(r) gets determined when the initial velocity 
profile v(r) of the cloud is given, in addition to the density profile. Since the center 
of the cloud is taken to be at rest in any spherically symmetric profile, the leading 
term in the Taylor expansion of v(r) is order r or higher, and it then follows from 
(3) that the expansion for f(r) begins with a term that is order r 2 or higher. So we 
expand f(r) as 

oo 

f(r) = Y,fnr n - (30) 

n=2 

and assume f 2 ^ 0. Although we wrote the expansion for F(r) in terms of the 
density p(r), we have written the expansion for f(r) directly, and not in terms of 
v(r). This is done for the sake of simplicity. Once the expansions for F(r) and 
f(r) are given, v(r) can always be deduced using Eqn. (3). Also, the assumption 
/2 7^ is a natural one, because if the cloud is not marginally bound, it will require 
a fine tuning of the velocity profile with the density profile to make f 2 zero. Once 
again we are interested in the limiting value of O, and from Eqn. (8) we find that 
it depends on the functions (3, p and G(—p) in addition to the already introduced 
function rj. These functions in turn depend on F and /, and using the expansion 
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for F(r) and f(r) from Eqns. (17) and (30), a tedious calculation leads to a power 
series expansion for 0: 

6(r) = -^1)72 [Gir + Q2r 2 + Qsr 3 + 0(r 4 )] . (31) 

The coefficients QuQ2 and Qs are rather complicated functions of / and F. For 
instance, Qi is given by 

Similarly, Q2 and Q3 can be expressed in terms of coefficients of F and /, though 
we do not give their explicit forms here. Go, which is the limiting value of 0, is 
fixed by choosing a in (31) such that 60 is non-zero, and finite. So we get 
©0 = Qu « = 5/3, if Q 1 ^ 0, 

©o = Q 2 , a = 7/3, ifQi = 0, Q 2 ^0, (33) 
©o = <?3, « = 3, ifQi=Q 2 = 0, Q 3 ^ 0. 
Only the coefficients upto 0(r 3 ) will be needed. This form for 60 should be con- 
trasted with the form it takes when / = 0, as given in Eqn. (23). There it is the 
density profile which determines 6 , while here the profile for / also enters the 
picture. 

The limiting value of the function H(X, r) in (7) is given by 

H(X,0) = X+^L. (34) 



The roots equation has the same form as in Eqn. (25), with Go given by (33). The 
nature of the singularity will be determined by the first non-vanishing coefficient in 
the expansion (31) for 0. If Qi ^ then from (25) and (33) it follows that 

Xl /2 = \Qu (35) 
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which should be compared with Eqn. (27), the corresponding solution when / = 0. 
The singularity will be naked whenever Qi > 0, and from its form given in (32) it is 
apparent that suitable choices of F and / will yield positive values for Qi, although 
the result cannot be expressed as neatly as when / = 0. In general, Qi will not 
be zero, since it depends on two independent quantities - the velocity profile and 
the density profile. Hence it is not quite realistic to consider the cases arising from 
setting Qi = 0, so long as fs and F\ are non-zero. However, if Fi and fs are both 
zero, then Qi will be zero, and we have to consider Q 2 . 

The expression for Q 2 , when Fi and fs are both zero, is analogous to that for 

Qi- 

* - - ( 2C ™ (t - i) t 1 + m) + ¥ - f ) ■ (36) 

In this case, a = 7/3 and from the roots equation we get X^ 2 = 3Q 2 /4:; the 
singularity will be naked whenever Q 2 takes positive values. When the metric and 
density functions are assumed to be even, Qi and also Qs to be given below both 
vanish and only Q 2 survives. This is the case considered by Newman in [5]. 

If /4 and F 2 are also zero, then we have to consider Q 3 , and in this case Q3 is 
given by 

* - - w) ( 3< w (I - i) i 1 + id + f - f ) • ^ 

This time we have a = 3 and the roots equation reduces to the same quartic as in 
Eqn. (29), the only difference being that now £ = —Q 3 /F^ 2 . The condition for 
the occurrence of the naked singularity is the same as stated following Eqn. (29) in 
terms of £. 

If the first three derivatives of the density are zero at the center, and in addition 
if /3 = j\ = fs = 0, then the first non- vanishing term in expression (31) for G will 
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be order r 4 , giving a > 3, so that the singularity will not be naked. This represents 
the transition to homogeneity in the case when / is non-zero. Not only should 
the first three derivatives of the density be zero, but the second, third and fourth 
terms in the expansion for / should also vanish, for the collapse to end like in the 
homogeneous case. 

C. The collapse of a cloud starting from rest 

This is an important sub-case of bound collapse, / < 0, as it may be argued that 
in this case the collapse is truly a consequence of gravity. From Eqn. (3) it follows 
that f(r) is now determined in terms of F(r): f(r) = —F(r)/r. The evaluation of 
G is now simpler, as we note in Eqn. (4) that G(—fr/F) = G(l) = ir/2. Hence 
to(r) = nr 3 / 2 /2y/~F, and using this in the defining relation 6 = t^y^X/r^ 1 we get 

= 4r i(a- io/2- ( 37 ) 

This is simpler than the form given for G in (8). Using arguments similar to those 
given in the previous section, it follows that if the qth derivative is the first non- 
vanishing derivative of the density at the center, then 

_ nqF q ri 
°° " " 4F r3(«-D/2' (38) 

where q = 1,2 or 3. This is essentially the same form as in (23). Also, the roots 
equation is the same as (25). Hence the nature of the singularity is precisely the 
same as in the / = case, except for a minor change in the coefficient of Go- The 
results of the marginally bound case apply. If q = 1 or 2, the singularity is naked. 
The case q = 2 is the one considered by Christodoulou [4]. If q = 3 then the roots 
equation is a quartic, same as in Eqn. (29) and the singularity will be naked or 

5 II 

covered depending on the value of the constant £ = 3ttFs/4F . If q > 4 the 
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collapse ends in a black-hole. It is interesting that the collapse of a cloud starting 
from rest has essentially the same features as a marginally bound cloud. 

IV. Global visibility and strength of the singularity 

So far we have been concerned with finding out conditions on the density and 
velocity profiles so that the collapse ends in a locally naked singularity. That is, the 
collapse is visible to an observer in the neighborhood of the singularity. However, 
it is still possible in this case that the null geodesies coming out of the naked 
singularity do not actually come out to a far away observer, but are covered by 
the event horizon and fall again in the singularity eventually. In such a case, the 
locally naked singularity could still be hidden in a black hole and only the strong 
cosmic censorship is violated but the weak form of censorship is intact. Of greater 
interest is the occurrence of a globally naked singularity, that is, the null geodesies 
emanating from the singularity escape the collapsing cloud entirely and reach an 
observer at infinity. It was shown in [9] that if a locally naked singularity occurs 
with a root Xo, it will be globally naked if and only if the condition 

7/(r)A(r) < aX . (39) 

is satisfied over the entire dust cloud. Using the definition of r\ and A and the 
Einstein equation (2) this can be written as 

p(r)r 3 " Q < aX . (40) 

Thus given a density profile which leads to a locally naked singularity, one has to 
check that the smallest positive real root satisfies the inequality (40) throughout 
the cloud, at the initial epoch. Clearly, one can always construct profiles which will 
satisfy this condition (as we demonstrate below with the help of two examples) and 
as a result, globally naked singularities would arise. 
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Let us take an example where / = 0, a = 7/3. This is the marginally bound case 
having first derivative of the density zero at the center, while the second derivative 
is non-zero. The condition (40) becomes p(r)r 2 / 3 < 7X /3. We take the density 
profile to be p(r) = po — P(2)^ 2 /2 with p( 2 ) positive. By evaluating the maximum 
of the left hand side in (40) and using the expression (28) for the root we find that 
the singularity will be globally naked if and only if the condition p\ < l.llp( 2 ) is 
satisfied. 

In our second example we take a = 3, without restricting /. Then (40) becomes 
p(r) < 3Xq, and for a density function which is decreasing outwards, this will be 
satisfied if and only if p < 3X . In the quartic (29) we take £ = —40, so that the 
quartic yields positive values for X and the solution X = 4F . Since po = 3F , 
the condition for global nakedness again is satisfied. 

We now examine the curvature strength of the naked singularity, which pro- 
vides an important test of its physical seriousness. For a detailed discussion of 
the criteria on curvature strength of a singularity and their applications to naked 
singular spacetimes, including the naked singularity in Tolman-Bondi models we 
refer the reader to [9] and [14]. Here we recall the result from [9] that the naked 
singularity is strong if a = rjo, and weak if a < tjq. In the present paper we have 
i]o = 3, and three possible values of a : 5/3, 7/3, 3. 

Consider first the marginally bound case. If the first derivative of the density at 
the center is non-zero, then a = 5/3, and the singularity is weak. If the first density 
derivative is zero, and the second one non-zero, then a = 7/3 and the singularity 
is again weak. If the first two derivatives are zero, but the third one is non-zero, 
then a = 3 and the singularity is strong. Thus restricting to even density functions 
leads to weak naked singularities (as was reported in [4] and [5]) and dropping the 
evenness assumption allows for strong curvature naked singularities. Precisely the 
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same results hold for a cloud starting from rest. Analogous results hold when / ^ 0. 
That is, if either Qi or Q 2 is non-zero, then a < 3 and the singularity is weak. If Qi 
and Q2 are both zero, and Q3 is non-zero, then a = 3 and the singularity is strong. 
The case a = 3 was discussed in [8] , though using a different scaling. It was shown in 
[9] that these results on strength hold for non-radial non-spacelike geodesies as well. 
Moreover, it is also possible that outgoing non-radial null or timelike curves could 
terminate at the singularity. Hence these results are a generalization of Newman's 
calculations on the strength along radial null geodesies [5]. 

The study of inhomogeneous dust is an important prelude to the study of other 
equations of state. As we mentioned before, departures from the Oppenheimer- 
Snyder work on collapse of homogeneous dust come in the form of introducing 
inhomogeneities as well as changing the equation of state. It is clear that introduc- 
tion of inhomogeneities significantly change the qualitative nature of the end state 
of collapse from a black hole to a naked singularity. As a next step, it will be of 
interest to examine the role of the equation of state. 
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